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Abstract 

Traditional smearing or blocking techniques serve well to increase the overlap of operators onto 
physical states but allow for links orientated only along lattice axes. Recent attempts to construct 
more general propagators have shown promise at resolving the higher spin states but still rely on 
iterative smearing. We present a new method of superlink construction which creates smeared 
links from (sparse) matrix multiplications, allowing for gluonic propagation in arbitrary directions. 
As an application and example, we compute the positive-parity, even-spin glueball spectrum up to 
spin 6 for pure gauge SU(2) at j3 = 6, L = 16, in D = 2+1 dimensions. 
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I. INTRODUCTION 



Operator construction in lattice gauge theory relies on the use of smeared and/or blocked 
links to improve the overlap onto physical states. (See [1] for a review of SU(N) gauge 
theory in 2+1 dimensions.) Smearing serves to reduce the ultraviolet fluctuations, allowing 
for a better projection onto the physically smooth wave-functions of glueball (and other) 
states, and blocking serves to keep the operator size physical as the lattice spacing decreases. 
Unfortunately, traditionally smeared links are only defined along axes directions, allowing 
construction of at most 2(D — l)-fold symmetric operators, with the rotation symmetries of 
the lattice timeslice. In D = 2+l, these operators couple to states with spin \J\ congruent 
modulo 4, ie a trial spin operator may contain contributions from spins 0, 4, 8, etc, and a 
trial spin 1 operator may contain contributions from spins 1, 3, 5, etc 

aa 

. To remedy the 

situation, links along either arbitrary or diagonal directions have been constructed 2j, IJ], 
allowing for operators with higher-order symmetry. These have been successful at resolving 
the 0/4 and 1/3 ambiguities p, |6| in the 2+1 dimensional pure gauge theory spectrum. 



The "path constructor method" of 



creates blocked links along the diagonal lattice 



directions for each blocking level, which are then used to define links between arbitrary sites 
i and j in a given timeslice. While effective, care must be taken of the details of the path 
construction for each unique arrangement of i and j. The "matrix method" of 2j, |4| creates 
the matrix of Green's functions between any and all sites % and j from the inversion of a 
matrix formed from the link variables of the timeslice, with a "hopping parameter" included 
to control the effective smearing level. While a single algorithm constructs all the paths 
at once, the inversion is painfully slow compared to iterative blocking and care must be 
taken to suppress unwanted torelon contributions. Note, however, that the correlation of 
rotated segment operators shown in J4] is much smoother for the matrix method, implying 
a better approximation of rotational invariance. Furthermore, the matrix method appears 
;o be directly applicable to fermion propagator calculations with some minor modification 
4], 3, Is]. For these reasons, further investigation of the matrix method is warranted. 

In this article, we explore a new method of smeared link construction based upon the 
matrix method but which alleviates its primary difficulties. Matrix inversion is replaced 
with matrix multiplications, and torelons are suppressed efficiently within the construction. 
Preliminary results for SU(2) in D = 2+l at (3=6 are promising. Section HT1 reviews the 
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technique for making operators of arbitrary spin J. Section II I II describes how to construct 
superlinks between arbitrary sites using matrix multiplications. In Subsection IIV Al we 
build traditional operators for Polyakov loops and for spins and 2 glueball operators. 
Construction of 8-fold operators is demonstrated in Subsection IIVB[ and "clock" operators 
with 12-fold symmetry are in Subsection IIV CI We highlight the difficulties with higher spin 
identification in Section [V] and Section |VT] presents results on an L = 16 lattice at (3 = 6. 
Some considerations and suggested improvements follow in Section IVII| and we conclude by 
summarizing. 

II. CONSTRUCTION OF OPERATORS WITH ARBITRARY SPIN ON A 
SQUARE LATTICE 

To resolve the modulo 4 ambiguity in the spin of conventional operators, one needs to 
construct operators of arbitrary spin from gauge- invariant loops at relative angles 9 = 2n/n 
other than the 7r/2 commonly available on the square lattice. These operators will suffer 
from their own modulo n spin ambiguities determined by the n-fold symmetry employed, 
but knowledge is power, and control over the value of n lets one select which sector of the 
spectrum one resolves. Loops which are 7r-symmetric will resolve the positive-parity, even- 
spin sector, and loops which are not 7r-symmetric will be able to resolve the full spectrum. 

Consider a timeslice of D = 2+1 in the continuum. Let C be the trace of the gauge field 
along an arbitrarily shaped closed loop of contour C aligned to an axis. Copies rotated 
through an angle 6 are denoted Cg. Then, to construct a gauge invariant operator Oj of 
arbitrary spin J, we take the weighted average over all loops Cg, where the weight is the 
complex phase e lje : 



(The shape and orientation of the loop may influence the phase cancellations.) On the lattice 
timeslice, our choice of angles is restricted to the n-fold angles of our chosen operator 



3=1 

and our loops Cg. will no longer be exact rotational copies of each other. Following [3], we 
normalize all these loops relative to their root-mean-square norm, in order to give the loops 




(1) 




(2) 
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approximately equal vacuum expectation values, ie 

% - C 6j x yf (((*,)) j/ity , (3) 

where (■ ■ - )j denotes the average over angles j. In practice, we find the smaller loops return 
roughly equal vevs after the above normalization, while the vevs of the larger loops display 
greater variance. Vacuum-subtracted correlation functions are defined by 

Cj{t) = (Okt)\Oj(0)) - (OfyiOj) , (4) 

and are normalized to Cj(0) = 1. For complex- valued operators we take the real part of 
Equation (j4]) as our correlation function. Similarly, we can compute the normalized cross- 
correlation for states of different J: 

{6Ut)\dM) 

Cjjit) - ' = , (5) 

^/(oi\Oj)(o\\Oj) 

where the tilde over O implies the vacuum-subtracted operator. The timeslice overlap 
between the operators is identified by the value of t in Equation (jSJ). We restrict ourselves in 
this article to using loops which are 7r-symmetric — if we were to use closed loops which were 
not symmetric under a rotation by 7r, we could recover both even and odd spins and both 
positive and negative parity, otherwise phase cancellations restrict us to the positive-parity, 
even-spin sector of the spectrum. 



III. SMEARED LINKS FROM MATRIX MULTIPLICATION 

To construct our arbitrarily shaped loops, we need to evaluate smeared links between 
arbitrary sites in the timeslice. We start by defining the link supermatrix A4, indexed by 
site, of dimension L ^ 1 x whose elements, themselves SU(N) matrices (giving a total 

matrix dimension of NL D ~ l x NL ^ 1 ), are constructed from the timeslice links Uij from 
sites i to sites j as follows: 

. if sites i and j are not nearest neighbors 
M^={ (6) 
Uij if sites i and j are nearest neighbors 



Note that M.ji = U},. Previously [2|, |4J, Green's functions were calculated by noting that 



hi 

1 



/C= ^— w l + aM + a 2 M 2 + ... (7) 

1 - aM 



for some small parameter a which behaves much like a "hopping parameter" for a scalar test 
particle sampling the gauge field. The matrix Ai k consists of all paths of length k — for even 
L, the odd powers of Ai are traceless while the even powers of Ai have a trace consisting 
of all paths of length k which terminate at their starting site, including those which double 
back on themselves. The diagonal elements of Ai 2 are "unity paths", 4 for each site, ie 
TrM 2 /ANL 0-1 = 1. The expectation value of the spatial plaquette (£/□) may thus be 
written as 



Note that this formula corrects that given previously [4j, |9j, as there are actually 7 paths 
which contribute unity per lattice direction and each plaquette is counted 8 times (4 corners 
by 2 orientations). Expectation values for larger closed loops may be found from the traces 
of higher even powers of Ai, eg the vev of the spatial 1x2 rectangle (U^-j) will relate to 
the trace of Ai 6 — these might prove useful if incorporating an improved action [10] into the 
calculation. 

Rather than computing the inverse matrix /C, which suffers from torelon contributions 
coming from the terms in the infinite series with k > L/2, we compute the sequence of 
matrices Ai k for k G {l,k max }. ("Torelon contribution" is actually a bit of a misnomer — 
the real problem is that for a distance k > L/2 between chosen sites i and j, Ai k j may 
contain contributions from both ways around the periodic lattice. With a suitable reduction 
in the unity path contributions (below) necessary for a path the short way around to be of 
distance k > L/2, these effects may be ameliorated.) For k max > L/2, torelon contributions 
are determined by the locations of sites i and j — we may propagate up to a total distance of 
L along a pure diagonal before contamination sets in. In practice we take k max = L to allow 
maximum propagation and avoid contamination by judicious choice of superlink orientation 
when constructing operators, or we may simply remove operators ah finito which appear to 
have been contaminated. Note that our parameter a has disappeared — we may reintroduce 
a smearing level parameter by applying a number p of smearing iterations, conveniently 
written as 



where mask(At p )x simply removes those elements of (Ai p ) 3 which do not correspond to 
nearest neighbors, before constructing our matrices Ai k ' p (at pre-smearing level p). Blocking 




(8) 



M p+ 



<= \mask(M p ) x (M p f - hM p ]/A 



(9) 
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algorithms may be implemented by a different choice of masking. We are now almost ready 
to use our superlinks to construct our gauge-invariant loops. 

In order to reduce the count of unity paths contributing to our JH k (dropping the index 
p), first we form 

M 2 = M 2 - diag(M 2 ) , (10) 

where we have removed the diagonal elements of A4 2 . Then we construct our sequence of 
matrices iteratively: 

M 2k+1 = M 2k x M , (11) 
M 2k+2 = M 2k xM 2 - diag(M 2k x M 2 ) , (12) 

removing the diagonal elements of the even powers as we go. As these matrices are sparse, 
performing the multiplications may be executed quite efficiently, and storage requirements 
are modest — the superlink construction requires but a fraction of the memory needed to 
compute the cross-correlations of the binned operators. We may verify our construction by 
checking that the spatial plaquette equals 

(U§ 4 ) = ( TrM^/NL - 1 - 12)/8 , (13) 

indicating a reduction in unity paths from 28 to 12. If we examine the superlinks along axes 
directions (Figured]), we see that to propagate a distance / between sites i and j along an 
axis, we should take -M'j 2 , and to propagate a distance I off-axis we use .ML-, where I is 
defined to be the number of links in the shortest path between sites % and j. 

The last step before we can construct our operators is to perform the reunitarization of 
the superlinks. The elements of Ai k , consisting of the sums of the possible paths, are not 
going to be elements of the gauge group. A preliminary calculation using unreunitarized 
superlinks projected well onto the ground states in the various sectors, but the excited states 
were impossible to recover cleanly. For the case at hand, SU(2), the reunitarization is simple: 
just divide each SU(2) element by its determinant. That simple step accounts for more than 
half the computational time currently consumed. For general SU(N), the reunitarization 
involves either an iterative procedure to maximize the trace of the projection or repetitive 



calculations of eigenvectors to unitarize the elementspjj, 112] — either way, a lengthy part of 
the calculation. Explorations of whether the matrix method for superlinks can return a 
naturally unitarized smeared link are underway. 



IV. CONSTRUCTION OF POSITIVE-PARITY OPERATORS 



A. Construction of Traditional Operators 

Operators for Polyakov loops may be constructed in D — 1 directions, at sizes I for which 
L/Z is an integer, simply by taking the product of superlinks along an axis p,: 

L/l-l 

®P,V = II ^i+mlfl, i+{m+l)l(i ' (^) 
m=0 

see Figure [2J "Box" operators of size Z centered at site i may be formed as 

O l D = Ml+ifa i+w x M x M\_i^ x M \_ Wi , (15) 

where ft, v are orthogonal directions, as in Figure [3J "Bar" operators, longer along the fi 
axis than the P, may be built using the clock points defined below as 

®B p, = M-X4 U Xmli x M.Xm\ u X2 t x M.X2 h XlOj 

( 16 ) 

X -MxiOi, Xm2, X -^Xm2 i ,X8 i x 7^X8;, , 

where Xmli is the on-axis midpoint between X2i and X4^, Xm2/ is between X2i and X4^, 
and similarly for O l B Figure HI Operators of spin and 2 are constructed from the bar 
operators as the normalized sum and difference of O l B ^ ^ and O l B ^ e , ie 

O l n ,j = Ol 0/2 = (O l B ^±O l B ^)/2. (17) 

Note that these operators are 7r-symmetric, 4-fold operators, hence the subscript 4. Here is 
the first example of when our construction algorithm returns trial loops of which some are 
highly redundant, sharing more than one superlink in common with loops of a different size. 
Inspection to remove redundant loops is best done graphically, giving easy identification of 
which loops are not redundant. States with spins \J\ congruent modulo 4 will have equivalent 
phases and hence the greatest mutual overlap, ie spins 0/4, 2/6, and 1/3 if we were using 
traditional L-shaped loops, will be mixed. The state with the lowest spin might not be the 
lightest state in a channel, thus care must be taken when assigning spin quantum numbers to 
the calculated states. From established results[l|, |4|, we expect the gap between the ground 
states for the + and 4 + , and between the 2 + and 6 + , to be large enough that the lower spin 
dominates the lightest few states. 



B. Construction of 8-fold Operators 



To achieve a degree of symmetry greater than the traditional 4-fold operators, we develop 
an algorithm to select points in a timeslice which are at relative angles of approximately 
27r/8 by placing an octagon on the lattice and labeling the points at size I as A\ through Hi 
clockwise, Figure [5j Using these we may construct 8-fold 7r-symmetric operators, eg 

O l 8 = M d Bu Al x M d Al> Fl x M d Fh El x M% t Bi , (18) 

with d = d(l) chosen appropriately, Figure [6l allowing us to resolve the spin 0/4 ambiguity. 
In practice, midpoints are also defined (as for the bar operators of the previous subsection) 
so that contamination is avoided. To build operators 0\ j with even spin J, we take the 
weighted linear combinations as above, Equation States with spins | J\ congruent modulo 
8 will have equivalent phases and hence the greatest mutual overlap, ie spins 0/8, 2/6, and 
4/12 will be mixed, but the overlap between the spin and spin 4 should be minimal. In 
practice, we find that the trial 4 + operator is more than capable of coupling to + states, 
and much care must be taken to distinguish a good 4 + from the other trial operators. The 
same applies to the 2/6 sector. As these operators are themselves "bar" operators, we may 
also form a pair of 2 + operators, O l 8 _ bar 2 , from appropriate pairings of the O l 8 . 

C. Construction of Clock Operators 

As before, we develop an algorithm to select points in a timeslice which are at relative 
angles of approximately 27r/12; ie, we place an imaginary clock face on the lattice centered 
at site X with increasing sizes and select those sites closest to where the hour marks would 
be, hence the appellation "clock" ' to operators made from this set of points, Figure [7J We 
label these points at size I as 1; through 12;, clockwise naturally. We define 6 rhomboid 
shaped loops at angles 9j = jvr/6 for j e {0, . . . , 5}, eg 

O l 0j = M d Xt Wl x M d Wu 12i x M d 2u 2i x M d h x , (19) 

as in Figure [HJ To build operators O l 12 j with even spin J, we take the weighted linear 
combinations as above. States with spins \J\ congruent modulo 12 will have equivalent 
phases and hence the greatest mutual overlap, ie spins 0/12, 2/10, 4/8, and 6/18 will be 
mixed. 
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V. IDENTIFICATION OF HIGHER SPIN OPERATORS 



The operators constructed above for the higher spins 4 and 6 are only exact in the 
continuum — on a finite lattice Hamiltonian eigenstates must transform as an irreducible 
representation of the square group C iv . Thus, our trial operators for both + and 4 + 
states will belong to the lattice irrep A±, and our trial 2 + and 6 + will belong to A%. Using 
superlinks smeared over some width to construct the loops, the phase cancellations intended 
to be angular may in fact be radial, resulting, eg. in a trial 4 + operator coupling to an 



excited + state. As has been known for some time [13( , "ascribing a definite spin to a single 
lattice state is therefore not possible without further dynamical information" coming from 
the behavior of the state as the continuum is approached. 

Meyer and Teper [J] present a couple of methods to identify accurately the spin of a 
calculated Hamiltonian eigenstate, and the interested reader is directed there. As the focus 
of this article is on the improved superlink construction and not on the precise identification 
of the higher spin operators (which requires the continuum limit to be taken), we perform 
a heuristic spin identification by combining the methods mentioned above. After the vari- 
ational procedure has orthogonalized the candidate states in channels J and J = J + 4, 
timeslice overlaps are computed via Equation (jSj) and the J states with the least overlap 
are identified as such. Consequently, the masses presented here for the higher spins 4 and 
6 should be viewed as preliminary until a continuum limit calculation using the superlink 
operators has been performed. 

VI. RESULTS 

A. Calculation Details 

To evaluate the performance of our superlink construction, we calculate the above opera- 
tors for pure gauge SU(2) on an L = 16 lattice at /3 = 6 with 15,000 measurements taken once 
every 10 compound sweeps, both with (p = 1) and without (p = 0) pre-smearing. Thermal 
updates are done via the Kennedy-Pendleton heat bath algorithm [14j augmented with a 
4:1 ratio of over-relaxation sweeps [15j and global gauge transformations every 19 sweeps. 
After selecting a set of operators based on their auto-correlation functions, cross-correlations 
are computed and used in a variational procedure to extract the lightest few states in each 



channel Timeslice overlaps are computed to check the performance of the variational 
procedure as well as the spin decomposition of the final states. 

Zero-momentum states are constructed by averaging operators over the timeslice. Ef- 
fective masses are extracted from the normalized correlation functions for t > a via the 
formula 

amf(t) = log ( ^4^) ■ (20) 



Cj(t) 

In principle, one should include a term representing correlations the other way around the 
periodic lattice, but in practice such effects are swamped by the statistical noise for lattices 
with extent L t great enough, especially for the heavier states. Errors are estimated using 
the jackknife procedure [7|. 

Standard practice is to observe the effective masses for evidence of a "mass plateau", 
which is taken to be the mass of the lightest state contributing to that correlation function 
over that region in t. While the variational procedure mentioned above does indeed produce 
operators which are in some sense orthogonal, there are certainly still contributions from 
two states (or more) for many of the correlation functions. In order to extract the most 



17] to the sum 



information possible from our correlation functions, we also perform a fit 
of two exponentials 

Cj(t) ~ A 1 exp(-m^) + A 2 exp(-m^) , (21) 

for a range of t e {0,4}, a notoriously hard problem. That difficulty shows up in a couple 
of ways: first, the 95% confidence limits on the parameters returned by the fitting routine 
are often quite large, and secondly, when unconstrained, some correlation functions return 
fitted parameters that are essentially nonsensical for the given application. The first we 
can circumvent by using the familiar jackknife procedure, the second by constraining the 
parameters to be physically significant positive quantities. To accurately calculate the binned 
values, care must be taken that the binned fits are sorted to match the order of states 
returned by fitting the ensemble correlation function. These two approaches to extracting 
the masses should be seen as complementary. 



B. Calculated Spectra 

In Tabled] we compare the vacuum expectation values for the lattice plaquette and spatial 
plaquette computed directly and the spatial plaquette computed via A4 4 and .M 4 . Our di- 
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rectly computed plaquettes are well within errors of the accepted value [lj], and our superlink 
plaquettes agree perfectly with the direct spatial plaquette. For the remaining operators, we 
perform the fits as above. We will examine the full table of fitted parameters and effective 
masses for the Polyakov loops and the + box operators to demonstrate and evaluate the 
effectiveness of the fitting procedure — for the remaining operators we will display only the 
best-fit masses and errors. Sometimes the correlation function is best fit by a single exponen- 
tial; that condition shows up as equal amplitudes of 1/2 for the two states. Otherwise, the 
state with the larger amplitude is taken as the dominant state in the correlation function. 
Note that the sum of the amplitudes was not constrained to unity — that resulted from the 
fitting procedure itself. 

Table [II] displays the best fit parameters and the effective masses in lattice units for 
the Polyakov loops after applying the variational procedure. Effective masses which were 
negative or imaginary have been zeroed out in the table. Jackknife errors are given in 
parentheses after the mass values. The apparent four-fold degeneracy is not unexpected, as 
we have included both orthogonal timeslice directions as well as two pre-smearing levels in 
the variational procedure. We highlight the ground state and also what appears to be an 
orthogonal excited state in boldface. The ground state Polyakov loop agrees well with the 
accepted value lj]. 

The + box operator's parameters and masses are shown in Table IIHI We note that now 
the degeneracy is two-fold, apparently a result of combining the two pre-smearing levels 
of operators in one variational procedure. Selected states are highlighted in boldface — 
preference is given to that state of a pair with either (or both) a single exponential or a 
lower mass, except for the heaviest, which was chosen for its lower error. The ground and 
first excited states compare favorably with previous results j^]. 

Satisfied that the fitting procedure is working reasonably well, for the remaining operators 
we present only the best fit parameters of the dominant state. The + and 2 + bar operators 
are in Table IIVI and the 2 + 8-bar operators of orthogonal orientation are in Table |V] These 
masses compare well with the values given by Teper and Meyer and are seen to be self- 
consistent. That the 2 + states are best fit by masses slightly higher than conventional seems 
to result from the conventional act of selecting the second effective mass as the best estimate, 
which is, after all, a highly arbitrary procedure. The + , 2 + , and 4 + best fit masses for the 
8-fold operators are shown in Table IVIL For the 4 + we have selected the state with the 
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least overlap with the + (see Table IVIip . which is in fair agreement with the value by 
Meyer, excluding the extremely heavy state, but we should keep an open mind as to the 
identification, noticing a degeneracy with an excited + in the same table. The + , 2 + , 4 + , 
and 6 + best fit masses for the 12-fold operators are shown in Table IVIIII The + displays 
and odd splitting of the ground state, and we highlight both values. The excited 2 + states 
are somewhat different in mass from the other operators' estimates. The highlighted ground 
states for the 4 + and 6 + are selected on the basis of their overlaps with the + and 2 + , 
respectively, Tables HXl and IXl (Overlaps between the 0/4 and 2/6 sectors are minimal.) 
The masses for the 4 + and 6 + compare well with those found by Meyer |12j ]. 



VII. CONSIDERATIONS AND IMPROVEMENTS 



The use of a pre-smearing level in this context serves more to remove ultraviolet fluctu- 
ations from the lattice configuration than to control how "wide" the superlinks are, and a 
better approach to controlling the effective width at various orientations should be devel- 
oped. Returning to the technique of removing the diagonal elements of the even powers of 
M. in Equation ffTTj) . the relative weightings of the more directional paths with unity path 
contributions might serve better to control how directional the superlinks are. As currently 
defined, the superlinks in Figure [D get effectively wider as their orientation approaches the 
diagonal: superlinks along an axis span the narrowest rectangle, while superlinks along a 
pure diagonal span the full square of links in-between. These exceptionally wide superlinks 
might explain why radial rather than angular phase cancellations seem to dominate the 
higher spin operators. Retaining the unity path contributions would give more weight to 
the more direct paths along the diagonal. One might also consider using superlinks of length 
1 + 2 rather than the minimum / used above — such superlinks would be proportionately wider 
and might be necessary as the lattice spacing decreases to keep the operators on a physical 
scale. 

The remaining modulo 4 spin mixing results from the mass eigenstates on a finite 2+1 
dimensional lattice properly belonging to irreducible representations of the square group, 
so our trial operators are going to contain contributions from multiple irreps. However, in 
the continuum limit, as the lattice spacing a approaches zero, our trial operators should 
couple to states of pure spin and mass. As computer power increases, simulations closer 
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to the continuum limit become possible, and the ability to create operators with arbitrary 
rotational symmetries will prove useful. 

Including both pre-smearing levels in a single variational procedure seems to have intro- 
duced a two-fold degeneracy for most of the operators, and in hindsight should probably have 
been avoided. Intermediate inspection of the correlation functions with and without pre- 
smearing indicated a better projection onto the ground states for the pre-smeared operators, 
so some level of pre-smearing should probably be retained in future calculations. 

A direct comparison of the performance of this technique with more traditional smearing 
techniques has yet to be done, as well as investigating how the performance scales with 
dimension D and gauge group N. Such comparisons are forthcoming. The negative parity 
and odd spin sectors need to be explored with appropriately shaped operators before the 
true utility of matrix superlinks can be determined. Relating these superlinks to those 
constructed previously 0, 4] by calculating the effective a of the current truncated expansion 
might shed light on both matrix propagator methods as well as on a new angle of approach 
to simulating staggered fermions on the lattice. 

Finally, we note that the superlink matrix Ai need not be confined to a single timeslice — 
one may form a super-superlink matrix M. D from all the links of the given lattice configu- 
ration. By dimensional reduction, the super-superlink matrix in D dimensions corresponds 
to the timeslice superlink matrix in D + 1 dimensions. Essentially, M. D holds all the infor- 
mation of the lattice configuration in a compact and convenient notation, and the traces of 
its even powers may be related to the ensemble plaquette average as well as terms used in 
various improved actions lOj. In light of Equation (J9j), familiar iterative algorithms may be 
expressed (and calculated) using efficient matrix notation, and with a suitable modification 
of the updating algorithm, a single matrix structure could be used for both updating the 
lattice configuration and calculating the operators. 



VIII. CONCLUSION 



The application of the matrix method of superlink construction has been demonstrated 
to agree with traditionally smeared operators and with previous methods of superlink con- 
struction. The primary difficulties with the matrix method have been alleviated by direct 
matrix multiplication rather than matrix inversion. The even spin, positive parity spectrum 
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for (3 — 6, SU(2) in D = 2+l dimensions displays excellent agreement with accepted values 
up to spin 6. The use of operators less symmetric than the ones presented here would re- 
solve the odd spin and negative parity sectors. Further exploration of the matrix method of 
superlink construction is warranted, both to resolve the complete spin and parity spectrum 
for SU(N) in higher dimensions and potentially to evaluate fermion determinants in SU(N) 
with quarks. 
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Plaquette Values 


(Uo) 


<^> 






0.8248964(11) 


0.8245756(15) 


0.8245756(15) 


0.8245756(15) 



TABLE I: Plaquette values computed directly and with superlinks. 









Polyakov Loops 








Best Fit Parameters 


Effective Masses 


Ampl 


Massl 


Amp2 


Mass2 


EffMl 


EffM2 


EffM3 


EffM4 


0.689 


0.994(017) 


0.311 


2.011(142) 


1.216(1) 


1.086(02) 


1.028(05) 


1.074(15) 


0.858 


1.077(014) 


0.142 


3.685(537) 


1.221(1) 


1.065(02) 


1.077(06) 


0.947(13) 


0.930 


2.588(026) 


0.070 


0.717(043) 


2.262(2) 


1.557(09) 


0.943(23) 


0.917(61) 


0.959 


2.564(143) 


0.041 


0.484(199) 


2.312(2) 


1.479(09) 


1.149(28) 


-0.000(35) 


0.890 


2.615(327) 


0.110 


0.804(172) 


2.169(2) 


1.425(07) 


1.248(26) 


0.217(37) 


0.971 


3.241(832) 


0.029 


0.546(2.016) 


2.904(4) 


1.580(18) 


0.723(42) 


0.469(71) 


0.953 


3.562(1.264) 


0.047 


0.713(3.222) 


2.994(4) 


1.406(18) 


0.888(39) 


0.406(65) 


0.500 


2.652(215) 


0.500 


2.652(538) 


2.669(3) 


1.653(15) 


2.405(162) 


-0.000(158) 


0.595 


4.803(1.270) 


0.405 


1.413(1.614) 


2.269(2) 


1.459(08) 


1.414(35) 


-0.000(43) 


0.925 


2.593(027) 


0.075 


0.750(041) 


2.257(2) 


1.560(09) 


0.995(25) 


0.729(55) 


0.812 


1.048(014) 


0.188 


3.271(694) 


1.232(1) 


1.068(02) 


1.079(06) 


0.941(14) 


0.683 


0.996(016) 


0.317 


2.039(123) 


1.226(1) 


1.089(02) 


1.030(05) 


1.070(15) 



TABLE II: Best fit parameters and effective masses in lattice units for Polyakov loops after appli- 
cation of the variational procedure. 
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+ Box Operators 


Best Fit Parameters 


Effective Masses 


Ampl 


Massl 


Amp2 


Mass2 


EffMl 


EffM2 


EfTM3 


EffM4 


0.961 


1.216(052) 


0.039 


5.415(721) 


1.257(1) 


1.201(02) 


1.306(009) 


1.502(039) 


0.793 


1.702(189) 


0.207 


2.991(1.844) 


1.864(1) 


1.752(08) 


1.716(042) 


-0.000(042) 


0.500 


2.100(068) 


0.500 


2.100(709) 


2.103(2) 


1.997(12) 


0.000(297) 


0.000(1.554) 


0.500 


2.279(074) 


0.500 


2.279(074) 


2.283(2) 


1.998(15) 


3.319(429) 


-0.000(416) 


0.500 


2.259(070) 


0.500 


2.259(070) 


2.579(3) 


2.476(34) 


0.000(078) 


0.000(137) 


0.551 


4.464(738) 


0.449 


2.065(321) 


2.761(4) 


2.136(26) 


5.951(1.983) 


-0.000(2.334) 


0.500 


2.634(107) 


0.500 


2.634(107) 


2.638(3) 


2.175(26) 


0.000(276) 


-0.000(282) 


0.869 


3.469(348) 


0.131 


1.486(1.166) 


2.871(4) 


2.017(27) 


1.588(130) 


0.634(278) 


0.500 


2.741(022) 


0.500 


2.741(022) 


2.740(3) 


2.796(50) 


0.000(096) 


-0.000(087) 


0.996 


2.732(164) 


0.004 


0.023(1.263) 


2.669(3) 


2.078(24) 


0.626(045) 


0.333(069) 


0.500 


2.586(284) 


0.500 


2.586(594) 


2.591(2) 


2.098(23) 


2.670(323) 


-0.000(368) 


0.500 


2.164(108) 


0.500 


2.164(1.127) 


2.170(2) 


1.939(12) 


0.000(307) 


0.000(1.421) 


0.500 


1.945(041) 


0.500 


1.945(213) 


2.112(2) 


1.911(11) 


2.306(112) 


-0.000(108) 


0.913 


1.241(023) 


0.087 


8.344(2.595) 


1.333(1) 


1.228(03) 


1.304(010) 


1.486(041) 



TABLE III: Best fit parameters and effective masses for + box operators after application of the 
variational procedure. 
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Bar Operators 


0+ 


2+ 


Amp Mass 
0.518 1.086(059) 
0.500 1.906(049) 
0.996 2.169(320) 
0.628 1.341(055) 
0.599 6.174(1.974) 
0.737 1.260(011) 


Amp Mass 
0.500 1.932(044) 
0.500 2.152(171) 
0.501 2.145(212) 
0.535 1.999(228) 
0.500 2.447(077) 
0.903 2.969(185) 



TABLE IV: Best fit parameters for + and 2 + bar operators after application of the variational 
procedure. 



2 + 8-Bar Operators 



Dirl 



Dir2 



Amp Mass 
0.500 2.045(029) 
0.926 2.553(388) 
0.595 2.236(314) 
0.528 2.151(293) 
0.500 3.705(066) 
0.847 2.896(275) 
0.500 2.527(104) 
0.500 2.283(147) 



Amp Mass 
0.998 2.038(122) 
0.936 2.576(514) 
0.500 2.581(132) 
0.500 2.960(059) 
0.500 3.252(097) 
0.500 2.781(086) 
0.500 2.627(266) 
0.500 2.243(218) 



TABLE V: Best fit parameters for 2 + 8-bar operators after application of the variational procedure. 
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8-Fold Operators 



0+ 



Amp Mass 
0.653 1.148(045) 
0.500 1.820(172) 
0.859 2.128(038) 
0.500 2.041(125) 
0.795 1.969(135) 
0.859 2.513(088) 
0.854 1.282(084) 
0.759 1.284(050) 



2+ 



Amp Mass 
0.500 2.030(019) 
0.967 2.491(410) 
0.500 2.574(083) 
0.500 2.801(122) 
0.726 4.197(716) 
0.999 2.842(273) 
0.500 2.164(096) 
0.992 2.364(320) 



4+ 



Amp Mass 

0.845 0.961(155) 

0.500 2.045(037) 

0.500 2.373(218) 

0.500 2.586(091) 

0.500 2.445(074) 

0.791 4.360(761) 

0.500 2.345(077) 

0.773 2.107(063) 



TABLE VI: Best fit parameters for + , 2 + , and 4 + 8- fold operators after application of the varia- 
tional procedure. 
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8-Fold Overlaps 



t=0 



t=l 



0+ 



4+ 



0.891 
0.377 
0.503 
0.406 
0.400 
0.467 
0.860 
0.559 



0.118 
■0.640 
0.326 
0.389 
•0.455 
•0.348 
•0.416 
•0.192 



0.238 
-0.288 

0.495 

0.506 
-0.538 
-0.421 
-0.361 

0.068 



■0.144 
0.344 
-0.313 
-0.601 
0.610 
0.282 
0.326 
0.033 



■0.299 
0.460 
■0.542 
■0.847 
0.625 
0.504 
0.495 
■0.049 



0.158 
-0.137 
0.229 
0.170 
■0.150 
-0.172 
-0.195 
0.050 



■0.100 
0.442 
■0.228 
■0.264 
0.267 
0.266 
0.301 
0.109 



0.719 
■0.309 
0.384 
0.344 
■0.299 
■0.366 
■0.706 
0.455 



4+ 



0.235 0.021 0.051 -0.025 -0.063 0.040 -0.021 0.192 
-0.067 -0.102 -0.035 0.048 0.062 -0.017 0.074 -0.055 
0.129 0.034 0.065 -0.036 -0.070 0.035 -0.025 0.106 
0.085 0.054 0.056 -0.055 -0.080 0.025 -0.039 0.071 
-0.088 -0.064 -0.058 0.055 0.073 -0.026 0.046 -0.071 
-0.117 -0.041 -0.057 0.036 0.066 -0.031 0.031 -0.096 
-0.208 -0.067 -0.059 0.045 0.082 -0.039 0.053 -0.170 
0.168 -0.029 0.029 -0.000 -0.027 0.027 0.016 0.139 



TABLE VII: Timeslice overlaps for the 4 + 8-fold operators. Rows correspond to the + operators, 
and columns correspond to the 4 + operators. 
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12-Fold Operators 



0+ 



Amp Mass 
0.518 1.108(049) 
0.871 1.365(031) 
0.500 2.159(084) 
0.528 1.854(250) 
0.585 3.440(566) 
0.558 1.403(041) 
0.564 6.204(1.214) 
0.583 1.424(035) 
0.905 1.210(018) 
1.000 1.864(100) 



2+ 



Amp Mass 
0.500 2.052(028) 
0.500 2.261(038) 
0.770 2.404(211) 
0.500 2.715(053) 
1.000 3.058(448) 
0.500 3.479(156) 
0.627 2.228(249) 
0.739 2.064(352) 
0.925 2.952(582) 
0.936 2.329(313) 



4+ 



Amp Mass 

0.500 1.483(074) 

0.500 2.528(098) 

0.500 2.645(087) 

0.500 2.658(112) 

0.500 3.020(046) 

0.710 5.169(848) 

0.500 3.284(070) 

0.500 2.984(136) 

0.500 3.140(173) 

0.500 2.178(065) 



6+ 



Amp Mass 

0.500 2.122(043) 

0.500 2.557(025) 

0.500 3.021(110) 

0.570 2.510(276) 

0.999 3.108(872) 

0.500 3.522(374) 

0.500 3.779(1.038) 

0.500 2.902(043) 

0.500 3.331(108) 

0.500 2.213(107) 



TABLE VIII: Best fit parameters for + , 2 + , 4 + , and 6 + 12-fold operators after application of the 
variational procedure. 
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12-Fold Overlaps 



t=0 



0+ 



t=l 



0+ 



4+ 



0.804 -0.063 0.000 0.084 0.001 0.013 -0.034 0.053 0.025 0.520 
-0.481 0.406 0.140 -0.043 0.097 -0.106 0.113 0.017 -0.019 -0.338 
0.338 0.517 0.332 0.187 0.306 -0.219 0.158 0.140 0.125 0.210 
-0.308 0.312 0.217 0.148 0.375 -0.312 0.128 0.068 0.107 -0.239 
-0.279 0.249 0.007 0.010 0.013 -0.176 0.013 0.061 -0.057 -0.260 
-0.435 -0.227 -0.214 -0.209 -0.347 0.241 -0.105 -0.105 -0.117 -0.258 
0.484 -0.075 -0.054 -0.043 -0.208 0.134 -0.123 0.005 -0.003 0.403 
0.602 0.298 0.258 0.153 0.242 -0.135 0.107 0.082 0.134 0.426 
0.747 -0.204 -0.054 0.056 -0.046 0.051 -0.069 0.014 0.026 0.503 
-0.047 -0.523 -0.221 -0.031 -0.164 0.147 -0.137 -0.073 -0.016 -0.004 



4+ 



0.211 -0.011 -0.002 0.026 -0.003 0.004 -0.006 0.017 0.004 0.135 
-0.143 0.057 0.016 -0.018 0.012 -0.015 0.015 -0.004 -0.008 -0.097 
0.068 0.050 0.028 0.021 0.023 -0.021 0.013 0.017 0.009 0.043 
-0.071 0.020 0.012 0.006 0.022 -0.019 0.008 0.003 0.006 -0.047 
-0.074 0.043 0.004 -0.006 0.005 -0.012 0.006 0.004 -0.009 -0.053 
-0.123 -0.004 -0.010 -0.027 -0.016 0.012 -0.001 -0.016 -0.010 -0.080 
0.101 0.010 0.005 0.008 -0.005 0.004 -0.000 0.008 -0.000 0.065 
0.138 0.024 0.019 0.024 0.015 -0.011 0.006 0.017 0.010 0.090 
0.199 -0.026 -0.007 0.023 -0.006 0.008 -0.009 0.013 0.005 0.129 
0.009 -0.071 -0.022 0.001 -0.016 0.019 -0.016 -0.010 0.006 0.012 



TABLE IX: Timeslice overlaps for the 4 + 12-fold operators. Rows correspond to the + operators, 
and columns correspond to the 4 + operators. 
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12-Fold Overlaps 



t=0 



2+ 



t=l 



6+ 



0.651 -0.078 0.002 -0.095 -0.158 0.156 0.074 0.204 0.203 0.476 
-0.650 0.463 -0.018 0.207 0.179 -0.220 -0.125 -0.426 -0.224 -0.390 
0.208 -0.565 0.116 -0.383 -0.355 -0.004 0.119 0.407 0.116 0.110 
-0.383 0.704 -0.122 0.444 0.141 -0.205 -0.065 -0.629 -0.146 -0.172 
0.147 -0.283 0.074 -0.371 -0.488 -0.237 0.084 0.214 0.132 0.127 
0.028 0.009 -0.062 0.068 0.179 0.439 -0.067 0.029 0.002 0.013 
0.554 -0.400 0.015 -0.182 -0.020 0.248 0.066 0.462 0.196 0.355 
-0.620 0.122 0.061 0.053 0.139 -0.137 -0.112 -0.174 -0.257 -0.457 
-0.108 0.488 -0.117 0.323 0.306 0.078 -0.111 -0.320 -0.052 -0.041 
0.561 -0.070 0.003 -0.132 -0.228 0.082 0.087 0.173 0.211 0.437 



6+ 



0.076 -0.009 0.002 -0.010 -0.020 0.020 
-0.071 0.041 -0.002 0.020 0.024 -0.023 
0.013 -0.043 0.008 -0.025 -0.021 0.003 
-0.035 0.050 -0.008 0.031 0.020 -0.016 
0.009 -0.017 0.006 -0.018 -0.023 -0.006 
0.006 0.002 -0.004 0.002 0.004 0.011 
0.060 -0.034 0.003 -0.018 -0.016 0.022 
-0.070 0.017 0.002 0.009 0.018 -0.019 
-0.003 0.037 -0.007 0.019 0.016 0.000 
0.064 -0.007 0.001 -0.011 -0.020 0.015 



0.009 0.024 0.024 0.054 
•0.012 -0.041 -0.024 -0.044 
0.009 0.030 0.007 0.006 
■0.011 -0.042 -0.014 -0.020 
0.006 0.014 0.006 0.008 
0.000 -0.000 0.003 0.003 
0.010 0.037 0.020 0.039 
•0.009 -0.025 -0.024 -0.048 
•0.007 -0.024 -0.002 0.001 
0.008 0.020 0.021 0.047 



TABLE X: Timeslice overlaps for the 6 + 12-fold operators. Rows correspond to the 2 + operators, 
and columns correspond to the 6 + operators. 
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FIG. 1: (Color online.) Superlinks of various widths and orientations. The superlink along the 
central line contains contributions from all the links covered by the surrounding shaded region. 



FIG. 2: (Color online.) Operators for Polyakov loops at sizes 1, 2, and 4. 
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FIG. 3: (Color online.) Box operators at sizes 1, 2, 3, 4, 5, 6, and 7. 



FIG. 4: (Color online.) Bar operators at sizes 2, 5, and 7. 
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FIG. 5: (Color online.) 8-fold points at sizes 1, 2, 3, and 4. 




FIG. 6: (Color online.) 8-fold operators at sizes 1, 2, 3, and 4. 
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FIG. 7: (Color online.) 12-fold clock points at sizes 1, 2, 3, 4, 5, 6, and 7. Note the redundancy in 
corner points for size 1 and between sizes 3 and 4. 
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